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Abstract 



An adiabatic approximation to the selfconsistent collective coordinate method is formulated in order to 
describe large amplitude collective motions in superconducting nuclei on the basis of the time-dependent 
Hartree-Fock-Bogoliubov equations of motion. The basic equations are presented in a local harmonic form 
which can be solved in a similar way as the quasiparticle RPA equations. The formalism guarantees the 
conservation of nucleon number expectation values. An extension to the multi-dimensional case is also 
discussed. 



I. INTRODUCTION 



Large amplitude collective motions (LACM), such as fissions, shape transitions, anharmonic vibrations and 
low energy heavy ion reactions, are often encountered in the studies of nuclear structures and reactions. To go 
beyond the phenomenological models assuming some macroscopic or collective degrees of freedom motivated 
by the experimental facts and intuitions, many attempts have been made to construct theories that are able 
to describe the LACM on a microscopic basis of the nuclear many-body Hamiltonian. Especially, theories 
based on the time-dependent Hartree-Fock (TDHF) approximation have been investigated extensively [p]-p^ . 
The TDHF is a general framework for describing low energy nuclear dynamics accompanying evolution of the 
nuclear selfconsistent mean- field 0,|8| . A LACM corresponds to a specific solution of the TDHF equation of 
motion. Since such a solution forms only a subset of the whole TDHF states (Slater determinants), it is often 
called a collective path, a collective subspace, or a collective submanifold. The collective coordinates are then 
a set of small number of variables that parameterize the collective subspace, and the collective Hamiltonian 
is a function governing the time evolution of the collective coordinates. One of the main purposes of the 
LACM theories is to provide a scheme to determine the collective subspace and the collective Hamiltonian on 
the basis of microscopic many-body Hamiltonian. Although the studies of LACM theories are the vast field 
of research with many recent developments to different directions, realistic applications to nuclear structure 
problems are rather limited. In this paper, we would like to propose a new practical method to calculate the 
collective subspace. 

The adiabatic approximation has been often utilized for formulating the theory of collective subspace. 
Indeed, some class of LACM, such as nuclear fissions, can be regarded as a slow motion, thus justifying 
the adiabatic approximation. The adiabatic TDHF (ATDHF) theory is one of the most well known 

adiabatic theories and has been applied in some cases to realistic description of heavy ion reactions |^ . The 
ATDHF theory, however, accompanied a problem of non- uniqueness of the solution [QH]. Later, efforts to 
settle the non-uniqueness problem were made from different view points. The works of Ref. |^ emphasize 
the importance of the canonical variable condition and the analyticity as a function of collective coordinate 
for finding a unique solution. The proposed procedure relying on the Taylor expansion method has not been 
applied to realistic calculations. Another work pointed out that the collective subspace can be uniquely 
determined by using the next order equation of the ATDHF theory. It is clarified also that the adiabatic 
collective path of LACM becomes the valley line of the potential function in the multi-dimensional space 
associated with the TDHF states 0-^. Further, the adiabatic collective path can be defined by equations 
for a local harmonic mode at each point of the collective path. These developments are summarized in a 
consistent way in the formalism of Ref. Note however that the adiabatic theory of Ref . |§ ] relys on a 
multi-dimensional classical phase space representation of the TDHF determinantal statesO,|l8|P A realistic 
application of this theory has not been done yet except for the one to a light nucleus |19|]. Furthermore a 
problem of particle number conservation arises when applied to superconducting nuclei | IG] . 

Theories without the adiabatic approximation have been also developed within the TDHF framework. The 
early works of this direction are called local harmonic approximations p^Jl^ . Later, a set of general equations 
that can determine the collective subspace and the collective Hamiltonian were found and formulated in a 
consistent form known as the selfconsistent collective coordinate method (SCC or SCCM) The theory 
is purely based on the TDHF with no further approximation. The method also provides a concrete and 
practical scheme to solve the basic equations using a power series expansion with respect to the boson-like 
variables defined as a linear combination of the collective coordinates and momenta. The pairing correlation 
in superconducting nuclei is easily incorporated within the SCCM by adopting the time-dependent Hartree- 
Fock-Bogoliubov (TDHFB) equation in place of the TDHF, and the conservation law of the particle number 



IJ. Thanks to these features the SCCM 
orations in medium and heavy nuclei [|l6|. 



is consistently introduced in the basic framework of the SCCM 
has been applied for many realistic descriptions of anharmonic vi 
However, the expansion method may not be suitable for the large amplitude motions of adiabatic nature, for 
which change of the nuclear mean-field is so large that the power series expansion of the collective coordinates 
may not be justified. 

In the present paper, we try to combine merits of two approaches mentioned above, i.e. the SCCM and 
the adiabatic thoery in order to formulate a theory that provides a consistent and practical method easily 
applicable to realistic descriptions of the adiabatic LACM in superconducting nuclei. We achieve this aim by 
introducing an adiabatic approximation to the general framework of SCCM. Here we treat superconducting 
nuclei since the pairing correlations play essential roles in many cases, like spontaneous fission, tunneling 
between superdeformed and normal deformed configurations, and coupling between coexisting states with 
different nuclear shape (shape coexistence phenomena) . Although the use of the superconducting mean field 
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requires us to respect the particle number conservation, the SCCM allows a simple and consistent treatment 
of the conservation law. We also avoid the non-uniqueness problem by utilizing the principles similar to that 
of Refs. Furthermore, we shall show that the equations of the adiabatic SCCM thus formulated can 

be transformed to another set of equations that have a similar structure as the local harmonic approach to 
the adiabatic theories Therefore, the present formalism is not only an extension of the SCCM, but also 
succeeds some aspects of the recent adiabatic theories such as Ref. |^ . 

In addition to the general formulation (Sect.||), we present a practical scheme to solve the basic equations 
given in the local harmonic form for general classes of the many-body nuclear Hamiltonian (Sect.|ll]| and 
Appendix). These equations are given in terms of the matrix elements of the many-body Hamiltonian 
expressed by the quasiparticle operators, thus enabling ones to develop a straightforward coding of a numerical 
program to solve the equations. In this way, we provide a complete procedure to extract the collective subspace 
and the collective Hamiltonian. We also discuss a possible prescription to extend the formalism to the cases 
of the multi-dimensional collective motions (Sect.^). Conclusions are outlined in Sect.^. 

II. BASIC EQUATIONS 
A. The sec method for superconducting nuclei 

In this subsection, we recapitulate the basic equations of the SCC method jl^ in a way suitable for treating 
the superconducting nuclei. 

We introduce the TDHFB approximation to describe LACM in superconducting many-fermion systems. 
Here the time-dependent many-body state vector is constrained to a generalized Slater determinant, 

which is chosen as a variational wave function. Time evolution of \(f>{t)) is then determined by the time- 
dependent variational principle 

s{<i>{t)\^§i-H\m = o, (1) 

where the variation is given by S \cf>{t)) — aj^a^j \(f>{t)) in terms of the quasiparticle operators {aj^,, Oq} which 
satisfy the vacuum condition Ua \ (t>{t)) = 0- 

We assume that the LACM can be described in terms of the collective variables, i.e. the collective coordinate 
and momentum {q,p} that are variables parameterizing the TDHFB state vector .0 The whole space of the 
TDHFB state vectors can be parameterized by M x (M — 1) variables (M being the number of the single 
particle states) as shown by the generalized Thouless theorem [|l^,|l^. A set of the TDHFB state vector 
\(j){q,p)) forms the collective subspace in which the LACM can be properly described. One of the main 
problems we concern is how to determine the collective subspace on the basis of the TDHFB equations of 
motion. At the same time, we need to determine the collective Hamiltonian TL{q,p) that governs the equation 
of motion for the collective variables {q,p}. This is a general purpose of theories of LACM. 

When we apply the LACM theories to nuclei in the superconducting phase, a special attention has to be 
paid to the particle number conservation. Since the TDHFB state vector is not an eigenstate of the particle 
number operator N , one would like to formulate the LACM theory so that the particle number expectation 
value is conserved during the course of collective motion. This is a problem which is specific to the TDHFB, 
but not to TDHF for which the state vector is a number eigenstate. 

It is well known ll^] that the expectation value of a conserved observable is kept constant during the 
time-evolution of \(l)(JYj governed by the TDHF(B) equations of motion. In the case of the pairing problem, 

the TDHFB state vector violates spontaneously the symmetry with respect to the gauge rotation e~'^'^'^ , but 
a rotational motion related to the gauge rotation (often called the pairing rotation) emerges automatically 
to restore the gauge symmetry. Therefore, the LACM of superconducting nuclei, described by the TDHFB 
theory, necessarily accompany the pairing rotation, for which we introduce the collective coordinate, (p, the 
gauge angle, and the conjugate collective momenta, N, which represents the particle number [p3. Thus, 



* We focus our discussion on a case of single collective coordinate. A multi-dimensional case is discussed in Sect. IV 
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we are obliged to consider a collective subspace that is parameterized by the set of four collective variables 
{q,p,^,N}^ 

Let us now present the basic equations of the SCCM that determine the collective subspace |(/)(g,p, (p, N)) 
and the collective Hamiltonian T-C{q,p, tp, N). As discussed above, the variable (p is introduced to represent 
the gauge angle. This requirement is easily satisfied |l5| if one uses the following parameterization 

\cj,{q,p,^,N))=e-^'^^\<j,{q,p,N)), (2) 

where N is the number operator of particles. Here \(j3{q,p, N)) represents an intrinsic state that rotates in 
the gauge space. 

Basic equations of the SCCM consists of a canonical variable condition and invariance principle of the 
time-dependent Schrodinger equation (TDHFB equation in our case). The canonical variable condition is, in 
general, given by 

{(j)iq,p,ip,N)\l-^\^{q,p,ip,N)) = 

{(j)iq,p,ip,N)\-^\^iq,p,ip,N)) = 
tdp 

d 

,p,Lp,N)\ I— \(j){q,p,p,N)) 

{(j>{q,p,ip,N)\-^\cP{q,p, ip, N)) 
idN 

for the collective subspace parameterized by two sets of coordinates {q, ip) and momenta (p, N). Although S 
is an arbitrary function of {q,p-, </5, iV}, we choose 5* = which is appropriate for the adiabatic approximation 
1^. Then the canonical variable condition can be rewritten as equations for the state \(j){q,p, N)), 

{(biq,p,N)\i-^\cP{q,p,N)) ^ p, (4a) 

{,p{q,p,N)\-^\^{q,p,N)) =0, (4b) 
{(t>iq,p,N)\N\^{q,p,N)) =N, (4c) 

{Hq,P,N)\-^\cP{q,p,N)) =0. (4d) 
idN 

The third equation requires that the collective variable N is identical to the expectation value of the number 
operator. In other words, the particle number expectation value does not depend on the collective variables 
{q,p) for the LACM under consideration. This is nothing but the condition of particle number conservation. 
The collective Hamiltonian is defined as value of the total energy on the collective subspace, given by 

n^Wq,P,ip,N)\H\^iq,p,ip,N)) (5a) 
= {<j){q,p,N)\H\^{q,p,N)). (5b) 

Since the Hamiltonian H commutes with the number operator N, the collective Hamiltonian does not depend 
on the gauge angle ip. Therefore, ip becomes cyclic as we expect. 

The invariance principle of the TDHFB equation plays a central role to determine the collective subspace, 
which requires that the TDHFB state vector \(j){q{t),p{t), (p{t), N{t))) evolving in time within the collective 
subspace should obey the full TDHFB equation, Eq.(|l|). This is equivalent to a condition that the collective 
subspace is an invariant subspace of the TDHFB equations of motion. Inserting Eq.(^) into the time- 
dependent variational principle, Eq.(p|), one obtains 



^ For simplicity, here we assume a single kind of particles. Extension to systems with many kinds (e.g., protons and 
neutrons in nuclei) is straightforward. 
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5 {^{q,p, N)\ H-f^P+f^Q + ^e-^N \cj^{q,p, N)) = 0, (6) 

where the infinitesimal generators defined by 

P\4>{q,p,N)) =i^|(/.(g,p,iV)), (7a) 

Q\<P{q,p,N)) ^\-l-^\<jy{q,p,N)), (7b) 

e|0(q,p,iV)) ^\^\<j,{q,p,N)) (7c) 

are used. These operators are one-body operators which can be written as hnear combinations of bihnear 
products {al^a^ij^apacnal^ap} of the quasiparticle operators defined with respect to \(j){q,p^N)). Because of 
the canonical variable conditions, these infinitesimal generators satisfy the following commutation relations 

(0(g,p,iV)|[Q,P]|0(g,p,iV)) =i, (8a) 
(0((?,p,7V)|[e,iV]|0(g,p,7V)) (8b) 

o o o „ 

and commutators of other combinations of Q, P, 0, N give zero expectation value. By taking the variation as 

o o o „ _ 

5 \4>(ci,p, N)) = {P, Q, 0, N} \(f>{q,p, N)), Eq.(P) produces the canonical equations of motion for the collective 
variables 

S " ^ " * (</'(9,P, N)\ [H, Q] \ct>{q,p, N)} , (9a) 
§ = =^mq,P,N)\ [H,°P] \cf>{q,p,N)) , (9b) 

J = 1^ = * {Hq,P,N)\ [H, e] \<l>{q,p,N)) , (9c) 

dt dtp ^ ' 

Using Eq.(|^), Eq.(^ then reduces to an equation of collective subspace 

6{cj,{q,p,N)\H-^P-^Q~-^N\cl,{q,p,N))=Q. (10) 

o o o ^ 

If we take a variation 5± that is orthogonal to the infinitesimal generators {P, Q, 8, N}, one can immediately 
show {(l){q,p, N) \ H \ (j>{q, p, N)) = 0, which implies that the the energy expectation value is stationary on 
the collective subspace for all the variations except for the tangent directions along the collective subspace. 
In other words, the collective mode is decoupled from the other modes of excitation. 

We remark here that the above basic equations of the SCCM arc invariant under point transformations of 
the collective coordinate 

q^q'^q'iq), (11a) 
p p' ~ p X {dq' /dq) \ (lib) 

The basic principles, i.e. the canonical variable condition, Eq.(^, and the invariance principle of the TD- 
HFB equation, Eq.(^), are not affected by the general canonical transformations of collective variables 
{q,p, ip, N} — > {q' ,p' , If' , N'}. By taking the parameterization, Eq.(^, and the specific choice of 5 = 
in Eq.(^, the allowed canonical transformations are restricted to the point transformations |^. 



B. Adiabatic approximation 

Assuming that the LACM described by the collective variables {q,p} is slow motion, we here introduce 
the adiabatic approximation to the SCCM. Namely we shall expand the basic equations with respect to 
the collective momentum p, which is appropriate for small value of momentum. Since the particle number 
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variable TV is a momentum variable in the present formulation, we also expand the basic equations with 
respect to n = N — Nq, when we consider a system with particle number iVo. 

Let us first consider the expansion of the TDHFB state vector \(j){q,p, N)) in the collective subspace. The 
origin of the expansion is the state \(t>{q)) = \(f>{q,p, N)) |p=o,JV=JVo- We can assume that this is a time- 
even state, i.e., T \(j){q)) — !</'('?)) under the time-reversal operation T (Here we consider the system of even 
numbers of particles). Thanks to the generalized Thouless theorem, the state vector \(j){q,p, N)) is expressed 
as 

\cl>{q,p,N))^e^'^'^'>'P''^^mq)) (12) 
by using the unitary transformation e*'^^'?^^'"). Here the hermitian operator G is given by 

G{q,p,n) = ^ (^Gap{q,p,n)ala}p + G*^i3{q,p,n)ai3a^ ^G{q,p,n)^. (13) 

Q>/3 

Here and hereafter, the quasiparticle operators {aJj,aQ,} are always defined locally at each value of q and 
satisfy the condition Cq, \4>{q)) = 0. We now expand the operator G{q,p,n) in powers of p and n and keep 
only the lowest order term. Namely, 

G{q,p,n)^pQ{q)+ne{q), (14a) 

Qil) = J2 {QM^UI + g;^(q)a/3a„) = Q(<z)^ (14b) 

Q>/3 

0(9) = E (0-/3(9)44 + 0a/3('?)«/3«") - 0(9)^- (14c) 

a>f3 

If we require that time-reversal of \(l){q,p, N)) causes sign inversion of the collective momentum p, i.e. 
T \(f>{q,p, N)) — \(f>{q,—p,N)), the operators Q{q) and 0(g) are time-even (TQ{q)T~^ = Q{q)) and time-odd 
{TQ{q)T^^ — — G)(q)), respectively. If we put n = (i.e. N = No), the parameterization Eq.(p^ together 
with Eq.(|lJ) reduces to \(j){q,p)) = e*^'^(9)j0(q)^ ^ which is the same form as the one introduced by Villars 
and often used in the ATDHF theories 
The collective Hamiltonian is expanded as 

Hiq,p, N) = V{q) + ^B{q)p^ + X{q)n, (15a) 
V{q) = n{q,p, N)\p^o,N=No = (0(9)1 H \^{q)) , (15b) 
Bil) = l^^^^^§^^\p=o,N=N„ = - im [[H,Q{q)],Qiq)] 10(g)) , (15c) 

m - ^^^|^|p=o.iv=Aro = (0(9)1 [H,^eiq)] 10(9)) , (15d) 

where we kept the collective momentum p up to the second order, while up to the first order in n. The 
collective Hamiltonian for the system with N = Nq particles ( n = ) is given by 

H{q,p,No) = V{q) + ^Biq)p^ (16) 

as the sum of the collective potential V{q) and the collective kinetic energy (the second term). 

We next expand the infinitesimal generators. It is convenient for this purpose to define the unitary trans- 

formation P' = e~*^Pe"^, Q' — e^^^Qe^^ ,Q' = e~*^6e*^ of the infinitesimal generators P, Q, 9. They are 
expanded as 

P' = P{q) + er^^^y^ = P[q) - - + (17) 



e 



-iG '-' iG 



e*^ - Q{q) + -[Q,pQ + ne] + (18) 
tap I 

= ^'^^4^/^ = ®(«) + ^[^'P^ + + (19) 

loN 2 
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with use of the general expansion formula 

g-iG^g.G ^ idG+ ^[idG,iG] + ^[[idG,iG],iG] + ... (20) 
The operator P{q) is the infinitesimal generator with respect to \(p{q)) defined by 

Piq)mq))^^§-^\Hq))■ (21) 

Similarly, we introduce the unitary transformation of the number operator and expand it as 

N' = e-'^Ne'^ ^ N + i[N,pQ + riQ] + ... (22) 
Substituting these operators in the canonical variable condition, Eq.(^), we have 

{^{q)\P'{q,p,N)\4>{q)) =p, (23a) 

{cl){q)\Q'{q,p,N)\cl,{q)) =0, (23b) 

{cPiq)\e'iq,p,N)\(l){q)) =0, (23c) 

{ct>{q)\N'{q,p,N)\4>{q)) ^N. (23d) 

Now we expand these equations with respect to momentum p and n. 
The zeroth order canonical variable conditions: 

(0(g)|P(g) 10(g)) = (0(g)|z^ 10(g)) =0, (24) 

(0(9)10(9)10(9)) =0, (25) 
(0(g)|e(g)|0((7))=O, (26) 
(0(g)| TV 10(g)) =7Vo. (27) 

Eqs-dH) and are automatically fulfilled by the definition, Eq. (|l^) , of the operators (5(g), 6(g). Eq.(|2^) 
can be satisfied if the g-dependent phase of |0(g)) is properly chosen. Eq.(p7|) is nothing but the constraint 
on 10(g)) for the conservation of average particle number. 
The first order canonical variable conditions: 

(0(9)1^ 10(9)) = -l> (28) 

(0(g)| [Q(g),e(g)] 10(g)) =0, (29) 

(0(g)| [Q(g), TV] 10(g)) =0. (30) 

One finds 

(0(g)| [Q(g),P(g)] 10(g)) (31) 



which can be derived by differentiating Eq.(|25D with respect to g and using Eq.(|2^). One can also derive 
from Eq.(|^ 

(0(g)| [P(g),iV] 10(g)) =0. (32) 

These equations give constraints on the infinitesimal generators Q{q), P{q) concerning the normalization, 
Eq.(|3l|), and the orthogonal condition to the particle number operator, Eq.(p2[). 

Next we expand the equation of collective subspace, Eq.(p^, to obtain a complete set of the basic equations 
for the adiabatic approximation. After rewriting Eq.(p^) as 

A - A dH. o dTi ° dTi ° 
S (0(g)| e-«iJe^« - ^P' - ^Q' - ^N' |0(g)) = 0, (33) 
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one can expand each term with respect to p and n with use of the equations hsted above. 
The zeroth order equation of collective subspace: 



S H - X{q)N - ^Q{q) \Hq)) = 0. (34) 

The first order equation of collective subspace: 

S {cj,iq)\ [H - X{q)N, Q{q)] - -B{q)P{q) |<^(g)) = 0. (35) 

I 

These equations are similar to the equations of path in the Villars' ATDHF theory except that the 
present paper deals with the superconducting Hartree-Fock-Bogoliubov (HFB) state, and that the Hamilto- 
nian accompanies the g-dependent chemical potential term —\{q)N. As we mentioned in Sect.|, the ATDHF 
theory has the problem that the solution satisfying these two equations is not uniquely determined . Al- 
though an additional validity condition was introduced to further constrain the solutions , the procedure 
of Ref. docs not fully solve the problem since the method does not work around the HF minima. 

The non- uniqueness problem has been investigated in recent developments of the adiabatic theories, and in 
our opinion they are classified into two different approaches. The first one represented by Ref. |^ claims that 
the solution is uniquely determined if an RPA boundary condition is specified at the HF minimum and if the 
analyticity of the collective path as a function of q is imposed together with the canonical variable condition. 
The solution, however, needs to be constructed in an analytic way or by means of a Taylor expansion method 
with respect to the collective coordinate q. We do not adopt this approach since we wish to construct a 
method applicable to systems under large excursion from the HFB minimum. We rather follow the other 
approach represented by Refs. 0-^. These theories require an additional condition that the equation of 
collective subspace (corresponding to the decoupling condition in Ref. ||]) should be satisfied up to the next 
order of the adiabatic expansion. In the present formulation, this second order condition is expressed as 
follows. 

The second order equation of collective subspace: 

5 {m\ \[[H - X{q)N, Qiq)],Q{q)] - B{q)AQ{q) \c^{q)) = 0, (36) 

where 

A<3(9) = ^+r(,)(5(<,), (37) 

This equation is equivalent in its mathematical form to the one given by Ref. if the chemical potential 
term —X{q)N is neglected. The last term —B{q)AQ{q), often called a curvature term, was simply neglected 
in the original version of local harmonic approximation p^, p^ . In the next subsection, instead of neglecting 
this curvature term, we shall rewrite AQ{q) and change Eq.(|36|) into a workable form. 

It is worth noting here the invariance of the adiabatic equations against the coordinate transformation. The 
collective momentum p undergoes the linear homogeneous transformation under the point transformation, 
Eq.(|l^). Therefore, different orders of the expansion with respect to the power of p are not mixed up under 
the point transformation. The invariance property of the basic equations of SCCM is thus inherited to each 
equation of the adiabatic approximation listed above. One can also confirm this property by seeing that the 
quantities appearing in the equations transform as 

Qiq) -> Q'{q') = Q{q{q')) ( ^] , (39a) 



P(g)^P'(g') = f'(g(g'))(^) \ (39b) 



dV dV _ dV 
dq dq' dq 



(39c) 
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(39d) 



AQ{q) ^ AQ'iq') = AQ{q). (39e) 

III. LOCAL HARMONIC APPROXIMATION TO COLLECTIVE SUBSPACE 

A. Local Harmonic Equations 

In this section we give an approximate but concrete procedure to construct a solution of the adiabatic 
sec method. To this end, we first derive, fi'om the adiabatic equations, another set of equations of collective 
subspace which can be solved in a way similar to the RPA equation. 

We first take a derivative of the zeroth order equation, Eq. (^) , with respect to q, which leads to 

S {<t>{q)\ [H - X{q)N, -P{q)] - C{q)Q{q) - ^^Q{q) " = 0, (40) 



^(9)-— -r(g) — , (41) 



d^V dV 

where AQ{q) and r(g) are given by Eqs.(^ and (|3^), respectively. Using Eqs.(^, we eliminate AQ{q) and 
rewrite Eq.(HG) as 

6 (0(g)| [H \iq)N, ip(<z)] - C{q)Q{q) ^[[^ - ^Qiq)],Qiq)] ~ |</.(g)) = 0. (42) 



Furthermore, due to Eq. (p4[) , we find 



^Q^iH~XN)A, (43) 



where {H — AiV)^ means a^a^ and aa part of the operator H — XN containing two-quasiparticle creation and 
annihilation in the normal-ordered expression. 

We thus replace Eqs.(|3^)-(|3^) by the equivalent set, 

S{<f>iq)\HM{q)mq)) ^0, (44) 
S {^{q)\ [HM{q), Q{q)] - \B{q)P{q) |0(g)) = 0, (45) 

5 {m\ [HM{q), -/{q)] C{q)Q{q) - ^^[[HM{q), {H A(g)iV)^], Q(g)] - |0(g)) = 0. (46) 

In Eq.dl) and (||), H - XN has been replaced by 

HM{q) = H-X{q)N-^Q{q), (47) 

since the last term has no influence. The operator HM{q) may be regarded as the Hamiltonian in the moving 
frame. The second and third terms can be identified with generalized cranking terms associated with the 
pairing rotation and the LACM, respectively. 

Equations (Es) and (Es) are linear equations with respect to the one-body operators Q{q) and P{q). They 



have essentially the same structure as the standard RPA equations except for the last two terms in Eq.(46). 
The quantity C{q) is the local stiffness parameter defined as the second (covariant) derivative of the collective 
potential V{q). The infinitesimal generators Q{q) and P{q) are thus closely related to the harmonic normal 
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modes locally defined for \4'{q)) and the moving frame Hamiltonian Hnjiq)- These equations may be called 
local harmonic equations. 

It was shown in Ref. that the zeroth, first and second order equations of ATDHF give a valley line 
of a potential energy surface in a multi-dimensional configuration space associated with the TDHF states. 
Similarly, the local harmonic equations we have obtained, Eqs.(H3)-(46), define the valley of the multi- 
dimensional potential energy surface. The solution of these equations will be uniquely determined if a 
suitable boundary condition is specified. These features are similar to the formulation of Ref. |^ where the 
valley equation of the potential energy surface is derived from the decoupling condition. 

We remark again that the local harmonic equations in the present paper differ from the ones of Rowe- 
Bassermann |ll2| and Marumori with respect to the the third and the fourth terms of Eq.(|46|), which 
arise from the curvature term (derivative of the generator) and the particle number constraint, respectively. 
It is important to keep the curvature term in order to maintain the relation between the collective subspace 
and the valley of the potential surface. We also note that the present formalism is invariant with respect to 
the point transformation of the collective coordinate, as is the formulation of Ref. S. 



B. Matrix Formulation of Local Harmonic Equations 

Let us now give a procedure to find the operators Q{q), Piq) that satisfy the local harmonic equations, 
( |45| ) and ( ^ ) , for a given state \c/){q)). Since they are linear equations with respect to these operators, it can 
be solved in an analogous way to the standard RPA. To show this, we first express the operator P{q) and N 
in terms of the quasiparticle operators: 

P{q) = t J2 (Pcp{q)aial ~ P:^{q)apa^ = P{q)\ (48) 



a>/3 



iV = 



(^"/3(9)«a4 + ^a/3(9)«/3«") ■ (49) 



Note that the a^a and c- number parts are neglected here since they do not change the state vector \(l){q)) except 
for the phase. The Hamiltonian H is also expressed in terms of the same quasiparticle operators. Assuming 
that the matrix elements Qaf3^ Paf3 are real, the local harmonic equations are written as the following matrix 
equations. 

{A - B)Q - B{q)P = 0, (50a) 

{A + B)P-C{q)Q- -^DQ-X'N = (50b) 
B{q) 

P^N = 0, (50c) 

2Q^P ^ 1, (50d) 

Here all quantities are functions of q, and Q = {....Qap, ■■■)'^ , P = {■■■,Paf:i, ■■■)'^ , and N = {...,Nap^ ■■■)^ , 
are the vector representation of the matrix elements with a > (3. A and B are the matrices whose elements 
are given by 

{A)a(),^S = 5a'i5ps{ea + Gp) + Wq^,^5, (51a) 

in terms of the matrix elements of the moving frame Hamiltonian 

Hniq) = ^eaajjOQ (52a) 

a 

+ i XI "a|,7'5"a'^/3«'5«7 (52b) 
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+ ^ 5Z («a!s7,5"a4"7"'5 + "A>/37"i"7«/3a«) • (52d) 



Here, due to Eq.(|3^), a^a^ and aa parts of HmW) vani sh, and a part of Hm{q) is diagonalized. The matrix 



elements of the residual interactions in Eqs.( 52b )-( p2cj ) are antisymmetrized with respect to the quasiparticle 
indices. The matrices A and B have the same structures as the ones often defined in the quasiparticle RPA 
formalism 

[0- The matr ix D is defined by 

(D)„;3,75 = I (0(9)1 mAliq), {H - \{q)N)A],aial + apa^^a^as] \^{q)) . (53) 
These matrix elements are expressed also in terms of the Hamiltonian matrix elements as 

(-D)a/3,7<5= {d^aP,-iS " d'^ap-fS + d]^^5l3S - dp\SaS " dl^Spj + dps6aj)/2, (54a) 

dfl3,jS^ Y{vl]3e,jhse - vD^.^gh^^ - vi^gh^e + vl^ghae), (54b) 
e 

d%^s^ Y^ivl^p^.M - vf^s^^h,^ + w^L,e^e/3 - vf^p.,K^), (54c) 

e 

C= II(""/375^75 - <L,/3/l75)> (54d) 
7>i5 

where hap is the matrix elements of {H — AiV)^ defined by 

{H - \N)A = '^hafi{a)a°^p + apaa). (55) 

a>/3 

Note that D contains the matrix elements of the type w^'^ and w'^^. These terms of the Hamiltonian do not 
contribute to the standard RPA equations. 

The solution of the matrix equations is obtained as follows. From Eq.(^), one obtains 

Q = \'B{q) {{A + B){A B) D- N, (56a) 

P = X'{A - B) {{A + B){A -B)-D- Viy^ N, (56b) 

with 

n - B{q)C{q). (57) 



The condition that the collective mode is orthogonal to the number operator, Eq.(50c), gives the following 
equation 

S{n) = N^{A - B) {{A + B){A B) -D- ^ly^ N = Q. (58) 

The quantity fl = B{q)C{q) represents the square of the frequency a; = ^/BC of the local harmonic mode, 
which is not necessarily positive. This equation can be regarde d as a dispersion equation to determine 
VI — up as a zero point of S{Vl). The normalization condition, Eq.( ^Ocj ), then gives a constraint on the value 
of X'^B{q). The value of the mass parameter B{q) is arbitrary, being related to the invariance under the 
point transformation Eq.(pl]). A choice of the coordinate, q, specifies a value of the mass parameter, B{q). 
In practice, the coordinate is often scaled so as to make the mass parameter unity. 

When the residual interactions are the separable forces such as the monopole pairing and the quadruple- 
quadrupole forces, the local harmonic equations reduce to a simpler form. The dispersion equation for the 
separable interaction does not require a matrix inversion as in Eq.([58[). The details are discussed in Appendix. 

Ref. has discussed a problem of spurious (Nambu-Goldstone) modes for local harmonic approaches, and 
stated that the RPA equation at non-equilibrium points must be extended in order to guarantee separation 
of the spurious modes. However, no practical way of solving the equation was given because the equation 
has parameters which we do not have a method to calculate. In our present formulation, the RPA equation 
is indeed extended to assure the number conservation. 
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C. Construction of Collective Subspace 



Let us finally give algorithms to construct the collective subspace \<p{q)) as a function of the collective 
coordinate q. Note that the local harmonic equations, Eqs.(^^-(^6|), are regarded as local equations in a 
sense that the equations can be solved independently for different values of q. At the HFB ground state, 
\(j>o), defined by the HFB equation 

6 (M H - XoN \M ^ 0, (59) 



we find dV/dq — 0. Therefore, l^o) is always a state on the collective subspace because Eq.(p4[) is automati- 
cally satisfied. Eqs. (^5|) and (^) reduce to the standard RPA equations at |0o) since the last two terms in 
Eq.(^) vanishes. The operators Q, P are then determined as one of the normal modes of the RPA equation. 
For non-equilibrium states, in general, Eq.Q) and the other two equation s, ji^ ) and (|6|), are coupled. We 



may solve the coupled equations in an iterative way. As discussed in Sect. [II B, one can find the operators 
Q{q)'^"'^ , P{q)^^^ by solving Eqs.(|45|) and ( ^6| ) for a given trial state 10(9))^"'' {n denoting the iteration step). 

This defines the moving frame Hamiltonian 77A/(g)^"+^' = H - X{qY"+^W - ^"^ Q(g)'"\ which can 

be used to construct a trial state 10(g)) for the next iteration. If the iteration converges, one obtains a 
state |(/'(q)) on which Eqs. (|^-(^ are simultaneously satisfied. Repeating the same procedure for different 
values of q, one finally obtains the collective subspace \4>{q)) and the collective Hamiltonian as a function of 

We remark here that the operator P{q) thus determined does not guarantee Eq.(pi|), although the other 
equations are satisfied. In this sense, the local harmonic solution is an approximate solution. The exact 



solution satisfying all the basic equations in Sect. 11 B may not exist in realistic situations. Only when the 



system is "exactly decoupled" g, the above procedure gives the exact solution. 

It is possible to choose another algorithm which satisfies Eq.(pT|) at the sacrifice of errors in Eq. (44). Let 



|(/)(go)) be a solution that satisfies the basic equations at g = gp. The infinitesimal generators Q{qo), Piqo) are 
determined by solving Eqs.(p5|) and (p6|). Then one can generate the state |0(go + Sq)) for an infinitesimal 
shift of the collective coordinate as 

|0(go+<5g)) =e-'^''^(*)|0(go)). (60) 

Repeating this procedure, one can construct a collective subspace. This solution should coincide with the 
one solved by the previous method if the system is exactly decoupled. Difference between the two gives a 
quality of decoupling for the collective subspace in the adiabatic approximation. The second procedure can 
be used to provide an initial guess, |(/)(g))^°\ for the iteration of the first method. 



IV. EXTENSION TO MULTI-DIMENSIONAL COLLECTIVE SUBSPACE 



In this section we extend the adiabatic SCC method to a case of multi-dimensional collective subspace 
described by D collective coordinates and conjugate momenta {q^,pi;i — 1, D}. 

One can easily derive the basic equations of the adiabatic SCC method in parallel to the derivation given 
in Sections by noting first that Eqs.(|l|Jl|) are now extended to 

|0(g,p,iV)) =e^^(«'P-") 10(g)), (61) 
G=p»Q'(g)+ne(g), (62) 

where the operator Q^{q) have now D components having coordinate label i. It is implied here and hereafter 
that the same coordinate index (i in the above expression) appearing in the super- and subscripts means 
to take the summation over it. The infinitesimal generator Pi{q) have also D components each of which 
is related to the derivative i-^ 10(g))- In the following, the coordinate dependence is often omitted. For 

instance, B^^{q) and Q^{q) will be simply denoted by i?*-' and Q*, respectively. 
The adiabatic collective Hamiltonian is expressed as 

n{q, p, N) = Viq) + ^B'^ {q)p,p, + \{q)n. (63) 
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The zeroth and the first order equations of the collective subspace are derived as 

6 (0(g) I H - X{q)N - I^Q' |0(q)) = 0, (64) 
6{ct>{q)\ [H^X{q)N,Q'] - ^B'^P, 10(g)) = 0, (65) 

while the second order equation becomes 

S Wq)\ ^[H- Xiq)N, Q\ + -^^Q\ Q\ Q'\ - ^J^'^'Q^ + ^''Q'-fe) I'/'C?)) = (66) 



with 



+ ^LQ^ (67) 



where B^j is the inverse matrix of B^^ and the bracket including three operators defined by 

= \{\\A,B\C\^\\A,C\,B\). (69) 

Expanding the canonical variable condition with respect to and n, the following equations are derived; 

(</.(g)|A|0(g)) = O, (70) 

(0(g)|iV 10(g)) = iVo, (71) 

and 

(0(g)| [Q',P,] 10(g)) =*<5,„ (72) 

(0(g)| [g\ TV] 10(g)) =0, (73) 

(0(g)|[A:,iV]|0(g)) =0. (74) 

These basic equations are invariant under the point transformation of the collective variables 

g' ^ g'^ = g'^(g), (75a) 
P^^P-=P, X (9gVag''). (75b) 

We have adopted the vector-tensor notation pO{ to manifest the transformation properties under the point 
transformation. Quantities which have a coordinate index in the subscript and in the superscript have the 
transformation properties of the covariant and contravariant vectors, respectively. For example, 

-> = X (ag'Vag^) , (76) 
P.^Pl = P, X (9gVag'0 . (77) 

The mass tensor _B*^ is the contravariant tensor of second rank. The operator Q^j defined by Eq.(p7|) is the 

covariant derivative of Q\ and F^^ is the ChristofFel symbol where the mass tensor Bij plays the role of 
metric tensor. 

Let us now derive local harmonic equations of collective subspace. Taking the g-derivative, the zeroth order 
equation (p3) leads to 



6 (0(g) I [H - X{q)N, -P,] - (g)Q^ - — Qf, - —N |0(g)) = 0, (78) 
d^V dV 
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As we have done for the D = 1 case, we would like to eliminate the covariant derivative Q^'^ in Eq.(|7q) in 
order to give a feasible form of the local harmonic equation. This was done for the D = 1 case with help 
of the second order equation of the collective subspace. The corresponding equations ( |66| ) give D{D + l)/2 
constraints, while number of unknown parameters, Q'^^, is D^. In fact, Eq.(|66|) is equivalent to 



.dV 



(80) 



where i?*-' are arbitrary one-body operators which are antisymmetric for exchange of indices i and j. If we 
choose R^^ — 0, we can eliminate the derivative term §^Q^i- Then, Eq.([78|) leads to 



S {^{q)\[H - X{q)N,^P^ - CjQ^ 



1 dX 

-[[H - X{q)N, [H - X{q)N)A].B,,Q^] - —^N = 0. (81) 




and is linear with respect to the i nfinit esim al gen erators Q^,Pi. We 



on the same line as in Sect. [II B and IIIC 



This equation is an analog of E- 
can numerically solve Eqs.(^ 

It should be remarked that the local harmonic equation Eq.(Bl|) for D > 1 are derived from equations, (|6j) 
and (^), but with an additional condition i?'-' = in Eq.(|8C). This condition is introduced to obtain the 
local harmonic equations parallel to the one-dimensional case. 



V. CONCLUSIONS 



We have formulated the adiabatic approximation to the general framework of the selfconsistent collective 
coordinate method in order to describe large amplitude collective motions in superconducting nuclei. The 
formalism, based on the TDHFB equations of motion, guarantees the conservation of particle number in 
a transparent way. We have shown that the equations of collective subspace are reduced to local linear 
equations for the infinitesimal generators, which can be solved with use of quasiparticle representation of 
the Hamiltonian matrix elements. This provides a complete procedure to determine the states e*?"^^?) \(j){q)) 
in the collective subspace and the collective Hamiltonian H{q,p) — V{q) + ^B{q)p^ as a function of the 
collective coordinate q and momentum p. A possible extension to the case of the multi-dimensional collective 
coordinates is also discussed. 

We emphasize that the equations given in this paper are solvable by means of the matrix method similar 
to the standard RPA. We hope that the present adiabatic theory is useful to solve number of open questions 
in the realistic studies of nuclear large amplitude collective motion. 



APPENDIX: SOLUTION FOR THE SEPARABLE INTERACTIONS 



In this appendix, we give solutions of the local harmonic equations of collective subspace for the case where 
the two-body interaction is given by the separable forces. We assume that the Hamiltonian is given by 

H = ho- ^F^F, (Al) 
where ha{— Kq) and F are one-body operators. Equivalently, one may write 

H ^ho~^F^+^F<^+W^F^-^F^-\ (A2) 

= (i^±F^)/2 ^i/'^^^t. (A3) 

For the separable forces, it is customary to neglect the Fock term of the forces. This approximation is easily 
and consistently implemented in the SCCM by assuming that the equation of motion for the time-dependent 
mean- field state \4>{t)) is now given by the time-dependent Hartree-Bogoliubov equation without the Fock 
terms, 

s{m^§i~Ht)\m^o, (A4) 

h{t) = ho^ (0(t)| F^+') \cb{t)) + kF^-^ (0(t)| F^-^ \cj){t)) . (A5) 
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The local harmonic equations (^4|)-(p6|) then become 

SWq)\hM{q)mq)) ^0, (A6) 

S [hM{q), Q{q)] - f^Q^F^--! - -B{q)Piq) |0(q)) = 0, (A7) 

S (0(9)1 [hM{q), ]B{q)P{q)] - - B{q)C{q)Q{q) ~ f^.+^'P^ 

~[F^-\ {h{q) - \{q)N)A]f'Q^ ~ InN \m) - 0, (A8) 
where Km (q) is the mean- field Hamiltonian in the moving frame defined by 

hAliq) = h{q) - ^Qil) - Kq)N, (A9) 

h{q) ^ho~ kF(+^ {^{q) I F(+) 10(g)) , (AlO) 
and definitions of other symbols are 

= Wq)\ [F(-), g(q)] |0(q)) , (Alia) 

= K {cb{q)\ [F^+\ \B{q)P{q)] \4>{q)) , (Allb) 

= -K {4>{q)\ [[F^+\ CHq) - X{q)N)A],Q{q)] |0(<z)) /2, (Allc) 

/iv = S(9)|^. (Alld) 

We express all operators in the above equations in terms of the quasiparticle operators {aj,, ad defined for 
hniq) and |0(<z)). For example, 

hniq) = ^eaaJ^Oa, (A12) 

F^^^ = E^S^(«^-4+ E^itl^«"«/'' (A13) 

Q!>/3 Q/3 

^^"^ = E^a?(«^-4-«^«-)+E^si4«/3. (A14) 

Q>/3 a/3 



We have assumed that all matrix elements are real. Equations ( A7 JA8) are then reduced to linear equations 



for the matrix elements Qai3,Pafj of the infinitesimal generators Q{q), P{q). They are easily solved to give 
the expression 

/ffi=tf'+/r. (A17) 

where we introduced the one-body operator 

i?(g)(±) ^ [Fl^\q), CHq) - A(g)iV)^] = ^ i?i^'(alat ^ a^a„), (AI8) 

Q>/3 



with Fl^\q) being the last terms of in Eqs.(|A13|jAl^. 
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Inserting this expression for the definition of /p^^ , /q \ we obtain the equations for unknown quantities 
f(+) A-) 



fpR^fg J In- Similarly, the condition of orthogonality to the number operator Eq.(50c) gives another 

equation for /p^\ /g 
a 3 X 3 matrix form 



equation for /p^\ /g \ /jv- They are summarized as a linear homogeneous equation which can be written in 



where 



Sxx' (^) 



oil — ^Op(+) p(+) -t- '^p(+) p(+) 



PR 

(-) 



0, 



Sl3 



(2) 

_F(+)F(- 



25 



(1) 

_F(+)_R(-) 



'fl(+)F(-) 



= 25 



(1) 



5 



521 = 25 



_F(+)Af 
(2) 

_F(-) F(+) ' 



(2) 



522 
523 



= 25 



(1) 

F(-) F(-) 



25 



(2) 

F(-)fl(-) 



1 



90(2) 
*-^F(-)Af' 



c _ o(l) 
"^31 — »^jvF( + ) ' 



532 
<S'33 



175 



(2) 

NF(-) 



(A19) 



(A20a) 

(A20b) 
(A20c) 
(A20d) 

(A20e) 

(A20f) 
(A20g) 
(A20h) 
(A20i) 



The functions S^^y with the symbols X,Y denoting {X,Y) = {F'^+\F^+^), {F^+\N), 

{N,N), (iV,F(+)), {N,R^-y) are given by 



5 



(1) 



E 

a>/3 



(Cq + 6/3)2 



f7 



0/3-1 Q/3, 



(A21) 



while the functions 5^|, with (A:,y) = , (i?(+',A^), (i^(-), F^), 

{F^-\N), iN,F'^-^) are given by 

4^1- = E 



c(>/3 



(ea + 6/3)^ 



SI 



:Ar„«Y^ 



The value of f2 is determined by finding the zero point of the dispersion equation 

det{5,,,(r!)} = 0. 



(A22) 



(A23) 



Normalizations of /p^\ /g \ /n are fixed by the condition Eq.(50d). It is straightforward to extend the above 
procedure to the case where the two-body interaction is given by a sum of the separable forces. 
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